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Abstract 

Supersymmetric unified models in which the Z' couples to the Higgs doublets, as 
in the Eq class of models, have large fine tuning dominated by the experimental 
mass limit on the Z'. To illustrate this we investigate the degree of fine tuning 
throughout the parameter space of the Constrained Exceptional Supersymmetric 
Standard Model (cEgSSM) that is consistent with a Higgs mass ~ 125 GeV. 
Fixing tan f3 = 10, and taking specific values of the mass of the Z' boson, with 
Mz> ~ 2 — 4 TeV. We find that the minimum fine tuning is set predominantly from 
the mass of Z' and varies from ~ 200 — 400 as we vary Mz> from ~ 2 — 4 TeV. 
However, this is significantly lower than the fine tuning in the Constrained Minimal 
Supersymmetric Standard Model (cMSSM), of 0(1000), arising from the large stop 
masses required to achieve the Higgs mass. 
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1 Introduction 



The Large Hadron Collider (LHC) has been accumulating data since 2009 with no obser- 
vation of new physics beyond the standard model (BSM) so far, placing strong limits on 
new coloured states in extensions of the standard model. For example, in supersymmetric 
(SUSY) models there are strong experimental limits on the first and second generation 
squark and gluino masses [TJ[2] which imply that they must be at least an order of mag- 
nitude larger than the electroweak (EW) scale. Within constrained versions of SUSY, 
where the stop masses are linked to first and second generation squarks masses, this can 
considerably increase fine tuning since the EW scale is very sensitive to stop masses, 
through the electroweak symmetry breaking conditions. 

At the same time Atlas and CMS have recently observed a new state consistent with 
a Standard-Model-like Higgs boson at m^ = 125 — 126 GeV [3j|4], which is within the 
range for it to be consistent with the lightest Higgs in supersymmetric models. In the 
minimal supersymmetric standard model (MSSM) this introduces further tension with 
naturalness since the light Higgs mass at tree-level is bounded from above by the Z 
boson mass (Mz)- The large radiative contributions from stops needed to raise it to the 
observed value typically imply very large fine tuning. For example the constrained MSSM 
(cMSSM) [5] has been shown to require fine tuning of (9(1000) if it is to contain a 125 
GeV Higgs mass (6j[7j . 

Here we consider fine tuning in an alternative class of constrained SUSY models which 
involves both an extra singlet field, denoted S, and an extra £7(1) gauge symmetry at 
low energy (TeV scale). As the singlet acquires a VEV, denoted s, it produces a ji term, 
denoted /x e ff, and it breaks the extra U{1) gauge symmetry, giving rise to a massive Z' 
boson. Such models can increase the tree-level physical Higgs boson mass above the 
Mz limit of the MSSM, due to both F-term contributions of the singlet and the D-term 
contributions associated with the Z', allowing lighter stop masses and hence reducing fine 
tuning due to stop loops. The exceptional supersymmetric standard model (E 6 SSM) [8j[9] 
is an example of such a model, inspired by the E$ group. At tree-level, the light Higgs 
mass is given as, 

m 2 , « M§cos 2 2/3+ — v 2 sin 2 2(3+— ^(1 + - cos 2/3) 2 + Am);, (1) 
"> ^ ' 2 4 4 
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where, tan/3 = — is the ratio between the two Higgs doublets' vacuum expectation values 
(VEVs), A is the Yukawa coupling of the singlet field to the Higgs doublets, and Am^ 
represents loop corrections. 

Indeed, Eq. [T] shows that the EgSSM allows larger tree-level Higgs masses than the 



NMSSM 10 , which in turn allows larger tree-level Higgs masses than the MSSM. This 
means that the E 6 SSM does not rely on such a large a contribution from the radiative 
correction term Ami m or der to reproduce the Higgs mass. As a result the E 6 SSM 
permits lower stop masses than either the NMSSM or the MSSM. In addition the A 
coupling in the EgSSM can be larger at low energies, while still remaining perturbative 
all the way up to the GUT scale, than is the case in the NMSSM. 
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One might conclude that this should lead to lower fine tuning in the E 6 SSM than 
either the NMSSM or MSSM, since the large stop masses are usually the main source of 
fine tuning in SUSY models. However, the origin of the extra term in Eq. [T] is due to D- 
terms arising from the coupling of the Higgs doublets to the extra U(l) gauge symmetry, 
and such D-terms also contribute to the minimisation conditions of the Higgs doublets. 
Indeed, as we shall discuss, one of the minimisation conditions of the E 6 SSM can be 
written in the form, 



M 2 Z _ „ 2 , K-m^ tan 2 /3) , ^Mf 
tan 2 /3 



-/4i+ v, i ;r": ^ +^ ; (2) 



where c, d are functions of tan which are of order ~ 0(1), m 2 d) m 2 are soft Higgs mass 
squared parameters and /x e g arises from the singlet VEV. Written in this form it is clear 
that the D-terms are a double edged sword since they also introduce a new source of 
tree-level fine tuning, due to the Z' mass squared term in Eq. [2j which will increase 
quadratically as M|,, eventually coming to dominate the fine tuning for large enough 
values of Mz>- This tree-level fine tuning can be compared to that due to yU e fr which 
typically requires this parameter to be not much more than 200 GeV, and similar limits 
also apply to Mz>- With the current CMS experimental mass limit for the Z' in the 



E 6 SSM of Mz> > 2.08 TeV 19 it is clear that there is already a significant, perhaps 
dominant, amount of fine tuning due to the Z' mass limit. 

In this paper we investigate this new and important source of fine tuning, namely 
that due to the Mz> limit, and compare it to the usual other sources of fine tuning in 



the framework of the Constrained EgSSM (cEgSSM) [13-16 . Although the impact of 
a SM-like Higgs with ~ 125 GeV on the parameters has recently been considered 
in [17, 18 , fine tuning was not considered. In fact the present study here is the first time 



that fine tuning has been considered in any supersymmetric E e model with a low energy 
Z '. To obtain the required Higgs mass in the cE 6 SSM, it turns out that the SM singlet 



field, S, must have a VEV s > 5 TeV as pointed out in 17 . This corresponds to a mass 



of the Z' boson predicted by the model of 1.9 TeV, which almost reaches the experimental 



bound of 2 TeV 19 . Thus, all the parameter space we study respects the experimental 
limit on Mz>- Fixing tan/3 = 10, and taking specific values of the mass of the Z' boson, 
Mz', ranging from 1.9 to 3.8 TeV we find that the current minimum fine tuning in the 
cEgSSM, consistent with a Higgs mass rrih ~ 125 GeV, varies from ~ 200 — 400, and is 
already dominated by the Mz> limit. However, this is significantly lower than the fine 
tuning in the cMSSM of (9(1000) arising from the large stop masses required to achieve 
the Higgs mass. 

The rest of the paper is organised as follows: Section two provides a short overview 
of the EgSSM. Then, the scalar Higgs potential and the electroweak symmetry breaking 
(EWSB) conditions are discussed in Section three. In Section four we discuss the fine 
tuning measure we use, and derive a fine tuning master formula for the E 6 SSM with a 
brief description of our semi-numerical procedure of calculating fine tuning. Section five 
is where we present our results and discussion, then we conclude the study in Section six. 
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2 The E 6 SSM 

The Exceptional Supersymmetric Standard Model (EgSSM) is a non-minimal supersym- 
metric extension of the SM, which provides a low energy alternative to the MSSM and 
NMSSM. It is well motivated both from more fundamental theories due to its connec- 



tion to E G GUTs, heterotic and F- string theory 12 and at the same time as a low 
energy effective model, providing solutions to phenomenological problems. For instance, 
as mentioned in the Introduction, the EgSSM allows a larger Higgs mass at tree-level 
than in both the MSSM and the NMSSM, thereby requiring smaller contributions from 
loops. In addition it also solves the \x problem associated with the MSSM by dynamically 
producing the /^-term at the TeV scale, without introducing the domain walls or tadpole 
problems that can appear in the NMSSM. 

The EgSSM is based on the Exceptional Lie group Eq. This contains both of SO (10) 
and SU(5) as subgroups, 

E 6 -»• SO(10) x U(1U (3) 
50(10) -> SU(5) x U(l) x , (4) 

and hence also contains the Standard Model gauge group, which is a subgroup SU(5). A 
linear combination of the two extra U(l)^ and U(l) x groups can survive to low energies, 
where it is spontaneously broken by a SM singlet field, S. This generates the mass of the 
associated Z' boson and the exotic quarks, as well as dynamically producing a fi e ff term. 
The model allows right-handed (RH) neutrinos to have Majorana masses at some scale 
between the GUT and low scales. This is achieved by choosing this linear combination 
to be, 



U(1)n = -^tf(l)* + \U(1) X (5) 

such that the RH neutrinos are not charged under U(1)n, hence it is possible to explain 
the tiny neutrino masses via seesaw mechanisms. 

At low energies, the group structure of the model is that of the SM, along with the 
additional U(1)n symmetry, 



E Q SU(5) x U(1) N (6) 
SU(5) ->■ SU(3) C x SU(2) W x U(1) Y (7) 

The matter content of the model is contained in the complete 27-dimensional represen- 
tation which decomposes under SU(5) x U(1)n to, 

27 — > (10, 1), + (5*, 2), + (5*, -3), + (5, -2) t + (1, 5), + (1, 0), (8) 

Ordinary Quarks and Leptons are contained in the representations: (10, 1) and (5*, 2). 
The Higgs doublets and exotic quarks are contained in (5*, —3) and (5, —2). The singlets 
are contained in (1,5), and finally the right handed neutrinos are included in (1,0). 

Moreover, the model requires three 27 representations, hence % = 1,2,3, in order 
to ensure anomaly cancellation. This means that there are three copies of each field 
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present in the model. However, only the third generation (by choice) of the two Higgs 
doublets, and the SM singlet acquire VEVs. The other two generations are called inert. 
Furthermore, in order to keep gauge coupling unification, non-Higgs fields that come from 
extra incomplete 27', 27 representations are added to the model. As a result, a jj term, 
which is not necessary related to the weak scale, is present in the model. 

The full superpotential consistent with the low energy gauge structure of the E 6 SSM 
contains includes both E e invariant invariant terms and Eq breaking terms, full details 
of which are given in M. However as in the MSSM it is necessary to forbid proton decay 
and therefore a generalsation of R-parity should be imposed, and additionally because 
the E 6 SSM includes three generations of every chiral superfield, there needs to be a 
suppression of new terms which can induce flavour changing neutral currents. To achieve 
this we impose either a Z^ symmetrj|^] (Model I) or a Z^ symmetry]^] (Model II) along 
with an approximate Z^ symmetry, under which all fields are odd except for the third 



generation Higgs superfields, which may arise from a family symmetry 20,21 
The Zif invariant superpotential then reads, 

+ l -M ij N c l N c 3 + fi'(H'W) + h%[E d H')e) + h%{H u H')N c 5 

+W MSSM ^ = 0), (9) 

where the indices a, (3 = 1,2 and i = 1, 2, 3 denote the generations. S is the SM singlet 
field, H u , and H d are the Higgs doublet fields corresponding to the up and down types. 
Exotic quarks and the additional non-Higgs fields are denoted by D and H' respectively. 

Finally to ensure that only third generation Higgs like fields get VEVs a certain 
hierarchy between the Yukawa couplings must exist. Defining A = A3, we impose K{ ~ 

^3 ^ ^1,2 3> fa/3, fa/3, hfj, h^j. 

3 The Higgs potential and the EWSB conditions 

The scalar Higgs potential is, 

V(H d ,H u , S) = A 2 |S| 2 (I^| 2 + \Hu?) + \ 2 \H d .Hu\ 2 

+ 9 -l(H ] d a a H d + Hla a H u )(Hla a H d + Ht<r a H u ) 

o 

+ 9 4(\H d f - \H u ff + 9 4(Qi\ H *\ 2 + Q^\ 2 + Qs\s\ 2 ) 2 (10) 

o Z 

+ m 2 s \S\ 2 + m 2 d \H d \ 2 + m 2 u \H u \ 2 
+ [XA x SH d .H u + c.c] + A Loops 

where, g 2} g'{= ^/3/5gi), and g[ are the gauge couplings of SU(2) L , U(l)y (GUT nor- 
malized), and the additional U(1)n, respectively. Q\ = — 3/-\/40, Q2 = — 2/v40, and 



4 All superfields except the leptons and survival Higgs are even. 

5 All the exotic quark, lepton and survival Higgs superfields are odd while all the other superfields 
remain even. 
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Q s = 5 /a/40 are effective U(1)n charges of H u ,H d and S, respectively. m s is the mass of 
the singlet field, and m u ^ = ^h u d - 

The Higgs field and the SM singlet acquire VEVs at the physical minimum of this 
potential, 



It is reasonable exploit the fact that s ^> v, which will help in simplifying our mas- 
ter formula for fine tuning as will be seen in Section 4. Then, from the minimisation 
conditions, 



9Ve & ssm _ 9Ve 6 ssm _ 9Ve 6 ssm 
dvi dv 2 ds 

the Electroweak Symmetry Breaking (EWSB) conditions are 



0, (12) 



Mf _ 1 (m 2 -m 2 tan 2 /3) (Q 1 - Q 2 tan 2 g) 

2" 2 AS+ tan 2 /3-l + 2 W + ^ + ^ s i tan 2 /3 - 1 

sin 2/3 « — ^ £ , (14) 



m 2 + ml + AV + ^-Q s s 2 (Qi + Q 



'2 



rnl™-\g'?Qy = -\Ml,, (15) 



where Mf = ifo' 2 + ^ 2 )( W 2 + v\) and Mf, « g[ 2 Q 2 s s 2 . 



Eq. 13 can be written in the form, 



Mf _ (m 2 -m 2 tan 2 /3) Mf, 

c-y - -A** + tan 2 ^ _ J + a-y. 

where c, d are functions of tan /3 which are of order ~ 0(1) and we have written [i c g = 
Written in this form it is clear that fine tuning will increase as Mz> increases. 
Another source of fine tuning is the large |/z e ff| term as mentioned in the introduction 



since satisfying Eq. 16 will require this term to compensate for any increase in either the 



second term (term 2: ~ m 2 ,m^) or the last term (term 3: ~ Mf,). 

The increasing experimental limits on M^/(~ s) results in constraining the parame- 
ter space of the E 6 SSM such that only relatively large values of mo and ma/2 result in 
successful solutions to the EWSB conditions (Fig. [j[ [IT]). 

Moreover, imposing universal boundary conditions, which is what characterises the 
cE 6 SSM, means that all low energy SUSY parameters can be expanded in terms of a few 
GUT-scale universal and fundamental input parameters, namely, 

m , ma/2, A Aj(0), /q(0), h tAr (0) (17) 

where, mo, ma/2 an d A are a universal scalar mass, a universal gaugino mass, and a 
universal trilinear coupling, respectively, and (0) means taking the parameter at the GUT 
scale (in the Results section, we refer to Aa(0) and ^,2,3(0) as Aq and Kq, respectively). 
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This is accomplished by using the one-loop RGEs of the scalar masses, so that one 
can express rn 2 Hu at the SUSY scale, Ms, as, 

m 2 Hu (M s ) = z x ml + z 2 m\ /2 + + z 4 m 1/2 A. (18) 
Then, it is possible to write, 

Ml A 2 

-^^j2 F ^ a i ( 19 ) 

i=l 

where, a denotes the fundamental parameters, z is the coefficient corresponding to each 
parameter, and is calculated numerically. F is some factor, possibly, involving tan/3. 

Whence, one can calculate (analytically or numerically) the sensitivity of Mz to each 
fundamental parameter, and this leads us to fine tuning. 



4 Fine tuning and the master formula 



To study the degree of fine tuning, a quantitative measure needs to be applied. Here 
we use the conventional fine tuning measure 22 , 23 , where the fractional change in the 



observable is calculated for a given fractional change in the input parameter, 



dh\M z 



d In a 



(20) 



where Mz is the mass of the Z bosorjj and a is one of the fundamental parameter in the 
set {m ,m 1 / 2 , A, A(0),k(0)}. 

For example, A a = 10 and 200 correspond to a 10% and 0.5% tuning in the parameter 
a, respectively. Moreover, for a given point in the parameter space, fine tuning is the 
maximum value of fine tuning in the set {A a }, and is denoted A max (or simply A). 

This measure has been used extensively within the literature e.g. (24 46 . 



4.1 Alternative tuning measures 

Some concerns have been raised in the literature regarding the use of this measure and its 
use is not universal with a number of alternative measures having been introduced and 
applied 47-61 . The {A a } measure the sensitivity of the parameters to the observable 



and as such are very dependent on the parametrization chosen. In particular whether one 
takes Pi to be the parameter or instead chooses a = pf introduces a factor two difference, 
and this factor two will then appear for every point in the parameter space. To remove 



this global sensitivity one can choose some normalisation 47-50 on the A a , however this 



then introduces questions about the bounds on the parameters and the probability is not 
clearly defined or understood. 

Additionally the overall tuning is chosen by taking A as the maximum of the individ- 
ual sensitivities {A a }, but a proposed alternative is to combine them in quadrature, like 



3 Note that some authors choose M§ instead of Mz- 



|A (Ml) 
in 0. 



Both measures can be easily linked since 
A a (M^). Our choice was made to enable straightforward comparisons with the results 
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uncorrelated errors 55 -58 . Clearly these measures can differ substantially, but it is not 
obvious which should be chosen. A new measure [60] defined tuning as the ratio of the 



parameter space volume (defined by fixed dimensionless variations in the parameters) to 
the same volume with the additional constraint that the dimensional variations of the 
observable are no greater than those of the parameters. As such this measure automat- 
ically combined the tuning from each parameter into a single tuning defined in terms of 
parameter space volume. For simple cases studied it was shown that this new measure 
was in greater agreement with the conventional measure than the alternative where the 
sensitivities are combined in quadrature, which might be understood as being due to large 
correlations between the individual sensitivities. 

Finally all the measures described so far define tuning as a theoretical feature of a point 
in parameter space, measuring how natural a point is. As such these measures quantify 
how natural phenomenologically acceptable points are once experimental limits have ruled 
out points which were initially favoured as being natural (or more natural). Instead 
within Bayesian analyses natural expectations for parameter space points, given by the 
prior distribution, are combined with experimental data to determine the probability 
defined as a degree of belief. If one must fine tune the parameters to get the measured 
values of observables correct, then this will correspond to only a tiny fraction of the 
total integrated prior volume, and therefore fine tuned scenarios should be automatically 
penalised. However in practice in MSSM studies Mz is often fixed to it's experimental 
value at the outset, reducing the dimensionality of the parameter space and missing the 
fine tuning. To fix this one can start off with a full set of parameters with the chosen 
prior distribution, unconstrained by EWSB requirements and then perform a Jacobian 
transformation |7, 62, 63j[65j|66]. The Jacobian factor accounts for the missed fine tuning 
and introduces similar derivatives as those appearing the sensitivity criterion, so it then 
appears as an effective "fine tuning prior". 

In the MSSM the conventional measure of fine tuning is numerically very close to this 
effective fine tuning prior (see e.g. |7]) and has sometimes been used directly as a fine 



tuning prior 33 67 , without directly calculating the Jacobian factor. 

Nonetheless the conventional tuning remains a very simple and useful measure and has 
continued to be used widely with the literature. We will employ it here for the following 
reasons: 

1. It is the most widely used tuning measure with which one can compare; 

2. It gives a good approximation of the effective fine tuning prior; 

3. It is simple to understand and apply; 



4. It provides a better match to the more complicated multi-parameter measure 60 
than combining sensitivities in quadrature. 

In particular please note that the simplicity and wide use is very important since this 
is the first quantitative investigation into tuning in this model and therefore comparison 
to what has been done in other models is of greater significance. Applying this measure 



7 This measure also allows one to combine several observables and had a normalised version of the 



tuning measure to deal with global sensitivity in a similar manner to 47 -50 , but with a slightly different 
normalisation and interpretation in terms of probabilities. 
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provides a quantification of the severity of tuning in the model, shows which regions 
have the least fine tuning and could be used as an "effective fine tuning prior" in future 
Bayesian studies of the model. 



4.2 Master Formula 



Having concluded the discussion on the motivation and suitability of this measure we 
now proceed to apply it in a quantitative analysis of fine tuning. To do so we first derive 
and present the master formula which gives the explicit expression from which the fine 
tuning is calculated. 



Using Equations [13 
for fine tuning in the EgSSM^J 



14, 15 and 20, we derive this master formula 



iC X 



9? 



Mf(tan 2 /3- 1) 



(1 - tan 2 0) «9(A V) + dm% _ ^ ^dm 2 u 



2 

tan (3 
cos 2/3 



da 

ds 2 AM 2 d 



da 



da 



da 



g 2 da 



(Q lC os 2 /3 + Q 2 sm 2 /3) 



1 + 



Ml 



m 2 + m 2 u + \ 2 s 2 + 9 -^Q s s 2 (Q 1 + Q 2 



(21) 



x 



X 



^ 9 ^ 1 - sin 2/3^(m 2 + m 2 + A 2 , 2 + |Wi + Q 2 )s 2 ) 



where 



1 - 



9? 



(tan 2 /3- 1) g 2 



(Qi - tan 2 PQ 2 ) x (Q 1 cos 2 f3 + Q 2 sin 2 0) 



(22) 



and f = (g' 2 + g 2 2 ). For tan/3 = 10; cT 1 ~ 0.88. 

The aim is to expand the low energy parameters, including s, in terms of the GUT- 
scale universal input parameters using the E 6 SSM RGEs as mentioned in the previous 
section. Next, the formula is implemented into a private cEgSSM spectrum generator 



(described in 15 16]) and fine tuning at each point in the scanned parameter space is 
calculated. In order to ensure accuracy of the results, the derivatives in the master 
formula for a — A(0) and a = k(0) are calculated numerically. And in order to calculate, 



da 



we use 



9?Q1 



(23) 



(24) 



where, as usual, m 2 is expanded in terms of the GUT parameters. 

Finally, throughout our study, we fix tan /3 = 10 since larger and smaller values restrict 
the availability of rrih ~ 125 GeV, and the parameter space 17 . 



Note we have left two terms in the second line of Eq. 21 written in terms of derivatives of cos (3 and 



sin j3 with respect to a. Substituting for soft masses here would unnecessarily clutter the expression and 
we note that these terms are numerically negligible since their contribution to fine tuning is very small 
(< 0(1)). This is due to the fact that they will be multiplied by an overall factor of order 0(< 10~ 12 ). 
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5 Results and discussion 



The scans are taken for fixed s = 5 — 10 TeV corresponding to Mz> = 1.9 — 3.8 TeV. We 
scan over 

-3 < A 3 (0) < and < /ti, 2 , 3 (0) < 3 

while fixing Ai^O) = 0.1 and tan (3 = 10. The other GUT parameters: mo,mi/2 and A 
are obtained as an output so that the EWSB conditions are satisfied to one-loop order. 
Then we plot both rrih and A max in the mo — plane. The key at the top-left of all 
plots corresponding to rrih shows the central value in a bin of width ±0.5 GeV, while that 
corresponding to A shows the central value in a bin of width ±50. 

Moreover, we select a benchmark point corresponding to each value of s. These points 
possess the smallest fine tuning in the mo — m\/2 plane consistent with a Higgs mass 
within the 124 < rrih < 127 GeV range, and m- g > 850 GeV. They are denoted as a black 
dot in Figures [[} [l2j These points and the relevant physical masses are summarised in 
Table [T] in Appendix [A| 




Figure 1: A max (left) and rrih (right) in the m,Q — m\ii plane for tan/3 = 10 and s = 5 TeV 
corresponding to Mz> = 1.9 TeV. We also fixed Ai^O) = 0.1 while scanning over —3 < A3(0) < 
and < ^1,2,3(0) < 3. The benchmark point corresponds to mo = 2020, = 1033 GeV. 



In the left panel of Fig. [T] the results for s = 5 TeV, corresponding to Mz> = 1.9 
TeV, are shown with fine tuning contours, ranging from 100 to above 800 for the highest 
mo- For each value of mo and m^, the parameters A, k, and A take different values. 
Since the Higgs mass strongly depends both on stop corrections and A, it will also take 
different values denoted by the Higgs mass contours displayed in the right panel of Fig. [TJ 
Since both fine tuning and the Higgs mass vary over the mo — TO1/2 plane the mass of 
the Higgs discovered at the LHC plays a crucial rule in fixing the level of tuning, though 
this dependence is significantly more complicated than in the MSSM. Thus, although for 
s = 5 TeV the tuning can in principle be as low as 100, in order to obtain ~ 124 
GeV the fine tuning must be more than twice as large as this. A benchmark representing 
points with the lowest tuning compatible with data shown as black dot in Fig. [T] having 
Abm = 251 with rrih ~ 124 GeV. Note that rrih ~ 125 GeV is almost impossible to 
achieve for s = 5 TeV (represented by the very small green region in the right panel). In 
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1.5 2 

m TeV 



1 1.5 2 

m TeV 



Figure 2: The left panel highlights the parameter responsible for the largest amount of fine tuning, 
^■max, in the m,Q — m.1/2 plane for tan/3 = 10 and s = 5 TeV corresponding to Mz' = 1.9 TeV. 
On the right a coarse scan shows which terms Eq. 16 give the largest contribution, with regions 
where the largest contribution comes from term 2, which is proportional to m\ — m\ tan 2 j3, are 
shown in yellow and while regions where the dominant contribution is from term 3, proportional 
to M|, are shown in blue. 



addition, the value M z < = 1.9 TeV slightly violates the CMS limit M z > > 2.08 TeV 19 



although this limit does not take into account the presence of lighter singlet states which 
increase the Z' width and reduce the leptonic branching ratio, weakening this limit as 



discussed in 13 



One also needs to take into account LHC constraints from squark and gluino searches 



which rule out m.1/2 J$ 1 TeV corresponding to a gluino mass m g < 850 GeV 17 . 

In Appendix [A] we provide a set on benchmark points corresponding to m.1/2 ~ 1 
TeV and these benchmark points are denoted by small black dots on the Figures. We 
emphasise that the cEgSSM has not been studied by any of the LHC experiments, and that 
the gluino mass limits in the E 6 SSM may differ from those of the MSSM as discussed 
recently |68|. Therefore, in choosing our minimum tuning benchmarks, the limits we 
assumed are quite conservative. From the results in [17], we find that in the cE 6 SSM the 
gluino mass is approximately given by m g ~ 0.85mi/2 and the first and second generation 
squark masses are given by trig ~ (1.3 — 1.8) mo, depending on m.1/2. In the future (for 
example when the full 8 TeV data set is analysed) the allowed values of m and m 1( / 2 
are expected to increase according to these approximate relations. Therefore, we show 
in Appendix [B] (Table [2]) the minimum allowed fine tuning associated with gluino mass 
in the 1 < m g < 1.5 TeV range, and the usual range for the singlet VEV s = 5 — 10 
TeV. Clearly, the fine tuning in the cEgSSM is not as large as that in the CMSSM, where 
increasing rrig to 1.5 TeV leads to minimum fine tuning > 1000 as found in [7j, while it 
varies between ~ 600 — 800 in the cE 6 SSM. 

At first sight, the distribution of fine tuning in the m — m 1( / 2 plane could seem counter 
intuitive since one might expect the region of smaller values of mo and mi/2 to possess 
lower fine tuning. However, the variation of A max can be understood by studying which 
parameter contributes the maximum fine tuning at each point in the parameter space. 
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We show this in Fig. [2] (left panel) where it is clear that the region of small mo and m 1( / 2 
is dominated by large fine tuning in the parameter Ao, resulting from a large |^ e ff| term 
in this region. 

In addition, k q can contribute to A max since A\ and m s are strongly dependent on 
this parameter. The physical origin of the fine tuning in k is due to the loops of exotic D- 
particles which serve to radiatively drive the singlet mass squared negative which triggers 
electroweak symmetry breaking. Finally, mo can be the source of fine tuning for very 
large values of mo which is the region extending beyond what we show in the plots. 

The relative fine tuning in the input parameters {m , m 1 / 2 , A, A(0), k(0)} does not 
directly tell us any information about the relative importance of the second and third 



terms on the right-hand side of Eq. 16 , both of which can independently be large and hence 
lead to a large \fi c s\ which is manifested as large fine tuning in A . It is therefore instructive 
to directly compare the magnitudes of the second and third terms of Eq. 16 where the 
former is proportional to vn? u and m 2 d , hence sfermions, and the latter is proportional to 
M§,. In Fig. [2] (right panel) we scan the parameter space for s = 5 TeV, and for each 
point we show which of the two terms is larger. The larger of the two would be responsible 
for the fine tuning at the corresponding point. It is clear, then, that Mz> (blue region) 
not only controls the minimum fine tuning allowed, but also is the dominating source of 
fine tuning over large regions of the parameter space. This is true for all the other values 
of s. However, some substantial contribution to fine tuning comes from sfermions as seen 
in the yellow region. 
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Figure 3: A max (left) and mh (right) in the uiq — mjrt plane for tan/3 = 
corresponding to Mz> = 2.3 TeV. The benchmark point corresponds to rriQ 
GeV. 



10 and s = 
- 1951, mi/2 



6 TeV 
= 1003 



As we increase s to 6 TeV (shown in Fig. |3J), we simultaneously satisfy the CMS 
mass limit on the Z' mass, with Mz> = 2.3 TeV, and we obtain more points with the 
heavier Higgs mass m^ = 125 GeV. Interestingly, the benchmark point in this case has 
a fine tuning Abm = 233 for w,h ~ 124 GeV which is slightly smaller than for the 
previous case with s = 5 TeV. Additionally, in the left panel in Fig. [3] a tiny region of 
A max = 200 appears as a small circle inside the A max = 300 band. While it is still Ao 
that is responsible for A max in that area as seen in the left panel in Fig. |4| this region is 
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Figure 4: The left panel highlights the parameter responsible for the largest amount of fine tuning, 
A m ax, in the m,Q — m^j plane for tan/3 = 10 and s = 6 TeV corresponding to Mz> = 2.3 TeV. 
On the right a coarse scan shows which terms Eq. 16 give the largest contribution, with regions 
where the largest contribution comes from term 2, which is proportional to m\ — m\ tan 2 /3, are 
shown in yellow and while regions where the dominant contribution is from term 3, proportional 
to M|, are shown in blue. 



associated with a slightly smaller \/i e e\ (|Ao|) and larger k q than in the adjacent regions, 
an effect which was not present in the results of s = 5 TeV. 

Moreover, Fig. [4] shows that the origin of fine tuning depends on the point in the 
mo — rrii/2 plane consistent with the Higgs mass and the LHC limits of squark and gluino 
masses, estimated above as m~ g ~ 0.85m!/ 2 and irtq ~ (1.3 — 1.8)m . For example if the 
squark and gluino masses are increased then it is possible that fine tuning is dominated 
by fine tuning in mi/2 or in Ao via large |/^ e ff| which could be due to heavy stop masses 
rather than large Mz' according to the right panel in Fig. |1} 



s = 7 TeV - Fine Tuning s = 7 TeV - Higgs Mass 




g, | 1 1 1 1 1 1 1 1 1 o | 1 1 1 1 1 1 1 1 1 
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m TeV m TeV 



Figure 5: A max (left) and (right) in the mo — m^ 2 plane for tan/3 = 10 and s = 7 TeV 
corresponding to Mz 1 = 2.6 TeV. The benchmark point corresponds to tjiq = 2186, mi n = 1004 
GeV. 
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s = 7 TeV - A max Parameter 
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Figure 6: The left panel highlights the parameter responsible for the largest amount of fine tuning, 
^■max, in the m,Q — m^j plane for tan /3 = 10 and s = 7 TeV corresponding to Mz> = 2.6 TeV. 



On the right a coarse scan shows which terms Eq. 16 give the largest contribution, with regions 



where the largest contribution comes from term 2, which is proportional to m d — m u tan j3, are 
shown in yellow and while regions where the dominant contribution is from term 3, proportional 
to M|, are shown in blue. 



For s = 7 TeV, corresponding to Mz> = 2.6 TeV, the region with ~ 125 GeV 
expands in comparison to s = 5 and 6 TeV, as can be seen by comparing the right 
panel in Fig. [5j to the previous plots. In addition a very small region with rrth ~ 126 
GeV appears for the first time. In the left panel of Fig. |5j fine tuning starts from 200, 
and reaches 600 outside the middle region. In addition, the tiny circle of points with 
smaller fine tuning than its surroundings in the small m — mi/2 region, which appeared 
previously in the results for s = 6 TeV, now grows a little. 

The chosen benchmark point has Abm = 270 for ~ 125 GeV. Notice how increas- 
ing s, hence Mz>, affects the lowest fine tuning possible in the parameter space, confirming 
that it is the Mz> term in Eq. 16 dominating fine tuning and defining its lowest value 
as can be seen in the right panel of Fig. |6j As before, this conclusion depends on the 
particular point in the m — "mxii plane. 
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s = 8 TeV - Fine Tuning 
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Figure 7: A max (left) and (right) in the uiq — m\j2 plane for tan/3 = 10 and s = 8 TeV 
corresponding to Mz> = 3.0 TeV. The benchmark point corresponds to rriQ = 2441, m-1/2 = 1002 
GeV. 
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Figure 8: The left panel highlights the parameter responsible for the largest amount of fine tuning, 
in the mo — m\/2 plane for tan /3 = 10 and s = 8 TeV corresponding to Mz> = 3.0 TeV. 



±max; 



where the largest contribution comes from term 2, which is proportional to mi — m^ tan z f3, are 



On the right a coarse scan shows which terms Eq. 16 give the largest contribution, with regions 

shown in yellow and while regions where the dominant contribution is from term 3, proportional 
to Mf , are shown in blue. 



For s = 8 TeV the Higgs mass ~ 125 GeV dominates over most of the itlq — ra\n 
plane as shown in the right panel of Fig. [7| Also the ~ 126 GeV region has become 
larger. However, fine tuning starts from 300, and the portion of the parameter space with 
A ma x > 500 is now more apparent than in the s = 7 TeV case. The Benchmark point 
has Abm — 302 for 125 GeV. The dominance of the Mz> term in Eq. 16 for fine 

tuning can be seen in the right panel of Fig. [8j with this conclusion dependent on the 
particular point in the m — m 1( / 2 plane. 
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s = 9 TeV - Fine Tuning s = 9 TeV - Higgs Mass 
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Figure 9: A max (left) and (right) in the uiq — mi/2 plane for tan/3 = 10 and s = 9 TeV 
corresponding to Mz' = 3.4 TeV. The benchmark point corresponds to too = 2709, m\/2 = 1001 
GeV. 
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m TeV m TeV 

Figure 10: The left panel highlights the parameter responsible for the largest amount of fine 
tuning, A max , in the mo — rrii/2 plane for tan/3 = 10 and s = 9 TeV corresponding to Mz> = 
3.4 TeV. On the right a coarse scan shows which terms Eq. [76] give the largest contribution, 
with regions where the largest contribution comes from term 2, which is proportional to m\ — 
m\ tan 2 (3, are shown in yellow and while regions where the dominant contribution is from term 
3, proportional to M|, are shown in blue. 



As we reach s = 9 TeV, corresponding to Mz> =3.4 TeV, which is shown in Fig. [9j 
we see that the region where ~ 125 GeV starts to shrink and is replaced by ~ 126 
GeV. If the Higgs mass is indeed ~ 126 GeV then there is a preference for s = 9 
TeV, especially for smaller values of tuq and m^. This illustrates the importance of an 
accurate determination in the Higgs mass for selecting the most appropriate value of s. 
Fine tuning starts from 200, although a very small region, and quickly increases to 500 
such that a significant portion of the parameter has A max > 500. The benchmark point 
has Abm = 330 for rrih ~ 125 GeV. The dominance of the Mz> term in Eq. [16] for fine 
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tuning can be seen in the right panel of Fig. 10, as usual dependent on the particular 
point in the mo — mi/2 plane. 




Figure 11: A max (left) and mh (right) in the mo — m-i/2 plane for tan/3 = 10 and s = 10 TeV 
corresponding to Mz> = 3.8 TeV. The benchmark point corresponds to mo = 2975,7711/2 = 1005 
GeV. 



s = 10 TeV - A m3X Parameter s = 10 TeV 
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Figure 12: The left panel highlights the parameter responsible for the largest amount of fine 
tuning, A max , in the mo — plane for tan/3 = 10 and s = 10 TeV corresponding to Mz> = 
3.8 TeV. On the right a coarse scan shows which terms Eq. [76] give the largest contribution, 
with regions where the largest contribution comes from term 2, which is proportional to m\ — 
m„ tan 2 @, are shown in yellow and while regions where the dominant contribution is from term 
3, proportional to M|, are shown in blue. 



Finally, for s = 10 TeV, corresponding to Mz> = 3.4 TeV, in the left panel of Fig. 11 



the fine tuning starts from 300, and the parameter space is severely restricted in terms 
of fine tuning as it is mostly covered by points with A max > 500. In addition, the region 
of rrih ~ 125 GeV has shrunk and now occupies a smaller portion than the rrth ~ 126 
GeV region. In addition a small region with ~ 127 GeV now exists prominently for 
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the first time (only a miniscule region existed for s = 9 TeV). Moreover, as seen before, 
the left panel in Fig. 11 contains short lines of points in the small mo — mi/j region with 
smaller fine tuning than their surrounding points for the same reason as before, namely 
that \fi e s\ can be somewhat smaller. 

The benchmark point has fine tuning Abm — 359 and m^ « 125 GeV. The dominance 
of the Mz> term in Eq. 16 for fine tuning can be seen in the right panel of Fig. 12, with 
the familiar dependence on the particular point in the mo — mi/2 plane. 



6 Conclusion 

Supersymmetric unified models in which the singlet VEV is responsible simultaneously 
both for /^ e ff and for the Z' mass, as in the Eq class of models for example, have relatively 
large fine tuning which is typically dominated by the experimental mass limit on the Z'. 
To illustrate this, we have investigated the degree of fine tuning throughout the parameter 
space of the cE 6 SSM. In fact this is the first time that fine tuning has been studied in 
any Eq model containing a TeV scale Z'. 

To quantify fine tuning we have derived a fine tuning master formula for the EgSSM 
and implemented it in a spectrum generator for the constrained version of the model. 
Using this we scanned the parameter space of the cEgSSM. The results are presented in 
the mo — mi/2 plane for fixed tan/3 = 10 and various s values corresponding to Mz> ~ 
2 — 4 TeV. This value of tan /3 = 10 is the optimum choice for achieving a large enough 
Higgs mass in the cE 6 SSM and so we have exclusively focussed on it here. We selected 
benchmark points corresponding to each value of s which possess the smallest fine tuning 
while allowing a Higgs mass within the 124 < m^ < 127 GeV range, and rrig > 850 GeV. 
They are the black dot points in Figures [!} [12} These benchmark points and the relevant 
physical masses are summarised in Table [T] for a gluino mass of about 900 GeV. Table [2] 
shows how the minimum fine tuning changes as the gluino mass limit increases up to 1.5 
TeV. As remarked earlier, the fine tuning in the cEgSSM is always significantly smaller 
than that in the cMSSM, for all gluino masses. 

It is clear that the Z' mass (determined by the s VEV value) has a significant effect 
on the naturalness of the cE 6 SSM model, with higher values leading to increased fine 
tuning. Therefore future improved direct mass limits on the Z' mass from the LHC will 
imply higher fine tuning. We have also seen an indirect relation between the Higgs boson 
mass and the Z' mass. For example if the Higgs mass turns out to be > 127 GeV 
then we are driven to s > 10 TeV corresponding to Mz> ^3.8 TeV requiring higher fine 
tuning. Conversely if the Higgs mass turns out to be m h < 124 GeV then s > 5 TeV 
corresponding to Mz> > 1-9 TeV allowing lower fine tuning. 

Given present limits, the results in Figures [!} 12 and Table [T] show that the present 



lowest value of fine tuning in the cE 6 SSM, consistent with a Higgs mass m^ ~ 125 GeV, 
varies from A ~ 200 — 400 where the allowed lowest fine tuning values, taking into account 
the relevant experimental bounds, are dominated by Mz> rather than the other sources 
of fine tuning. This is presently significantly lower than the fine tuning in the cMSSM of 
A ~ 1000 arising from the large stop masses required to achieve the Higgs mass. 

In the future, the LHC lower limits on gluino and squark masses will improve, along 
with the Z' mass limit (or else a discovery will be made) and the Higgs boson mass 
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will be more accurately specified. It is not completely clear where the dominant source 
of fine tuning in the cEgSSM will originate from in future. However the results in this 
paper allow this question to be addressed. The future Z' mass limit will determine the 
minimum s value permitted, while the Higgs mass and gluino and squark mass limits will 
determine the allowed regions of the mo — mu2 plane, from which the fine tuning may 
be read off from the contour plots we provide. 
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A cEeSSM Benchmark points 

Table [T] lists the details on the masses and parameters associated with each benchmark 
(BM) point that was chosen. We can see that m increases significantly as s {Mz>) 
becomes larger, while TO1/2 is roughly constant. Upon choosing a BM point, we imposed 
the limit m\ii > 1 TeV to have gluino mass m g > 850 GeV. The gluino masses for our 
benchmark points are about 900 GeV or close to it, hence if the experimental limits on 
■trig are to be increased for constrained models, then fine tuning will increase as well. The 
lightest stop, tx, masses range from 1.7 TeV to 2.4 TeV for the range of s we studied, and 
thereby is above the experimental limits. 
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BMl 


BM2 


BM3 


BM4 


BM5 


BM6 


s [TeV] 

L J 


5 


6 


7 


8 


9 


10 


tan j3 


10 


10 


10 


10 


10 


10 


X 3 (M X ) 


-0.228392 


-0.264564 


-0.25 


-0.237651 


-0.22603 


-0.217085 


Ai, 2 (M x ) 


0.1 


0.1 


0.1 


0.1 


0.1 


0.1 


Kl,2,3(M x ) 


0.17598 


0.192349 


0.211111 


0.228859 


0.245226 


0.260101 


mi/2 [GeV] 


1033 


1003 


1004 


1002 


1001 


1005 


m [GeV] 


2020 


1951 


2186 


2441 


2709 


2975 


A [GeV] 


-83 


500 


661 


781 


846 


888 


171^(1,2,3) [GeV] 


2252 


2234 


2659 


3149 


3680 


4222 


1714,(1, 2, 3) [GeV] 


3186 


3501 


3991 


4499 


5017 


5540 


^(1,2,3) [GeV] 


1782 


2238 


2752 


3279 


3812 


4347 


\m x o\ [GeV] 


1973 


2349 


2727 


3105 


3483 


3861 


m h3 ~ [GeV] 


1889 


2267 


2645 


3023 


3401 


3779 


m x p [GeV] 


1809 


2189 


2566 


2944 


3322 


3699 


m s (l,2) [GeV] 


2448 


2548 


2897 


3263 


3639 


4014 


m ff2 (l,2) [GeV] 


1970 


1847 


2023 


2218 


2426.5 


2633 


m Hl (l,2) [GeV] 


1887 


1685 


1824 


1986 


2167 


2343 


/J#(l,2) [GeV] 


492 


569 


642 


711 


777 


841 


m,a(l,2) [GeV] 


2505 


2461 


2687 


2934 


3199 


3468 


ra u - ~ rrir (1.2) [GeV] 


2553 


2507 


2729 


2973 


3235 


3501 


m d - 2 (l,2) [GeV] 


2571 


2558 


2810 


3082 


3372 


3665 


m el (l,2,3) [GeV] 


2136 


2107 


2366 


2641 


2935 


3224 


m e ~ 2 (l,2,3) [GeV] 


2267 


2271 


2550 


2848 


3159 


3468 


m fl [GeV] 


2119 


2090 


2347 


2623 


2912 


3200 


m^ 2 [GeV] 


2259 


2263 


2541 


2838 


3148 


3457 


m 6 [GeV] 


2202 


2151 


2340 


2549 


2777 


3009 


m 6 ~ [GeV] 


2552 


2539 


2789 


3059 


3347 


3639 


m t ~ [GeV] 


1741 


1681 


1839 


2016 


2212 


2411 


m t ~ [GeV] 


2215 


2166 


2354 


2561 


2787 


3018 


m v o ~ |m,±| [GeV] 


887 


1174 


1258 


1329 


1386 


1443 


mft 2 ~ ~ m#± [GeV] 


1890 


2268 


2646 


3025 


3403 


3782 


m h [GeV] 


124 


124 


125 


125 


125 


125 


mj [GeV] 


901 


879 


887 


892 


898 


906 


|m^±| ~ |m v o| [GeV] 

Ai 1 A2 1 L J 


285 


279 


279 


279 


279 


280 


m v o [GeVl 


162 


157 


158 


158 


158 


158 


^max 


251 


233 


270 


302 


330 


359 



Table 1: Parameters and masses for the benchmarks with lowest fine tuning and Higgs masses 
in the range of = 124 — 125 GeV in the cEqSSM. 



B Fine tuning and rrig 

As the lower limits on the gluino mass are expected to rise, Table [2] shows the minimum 
amount of the fine tuning corresponding to different values of gluino mass within rrig = 
1 — 1.5 TeV, and for s = 5 — 10 TeV. The corresponding Higgs mass is shown in parenthesis 
next to each value of fine tuning. 
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s [ievj 


5 


6 


7 8 


9 


10 


m~ [To VI 

nig [lev j 


A (m„ [GeV]) 


1 


293 (124) 


297 (124) 


324 (125) 


367 (125) 


405 (126) 


443 (126) 


1.1 


388 (125) 


348 (124) 


358 (124) 


408 (125) 


454 (126) 


497 (126) 


1.2 


474 (124) 


440 (125) 


400 (124) 


448 (125) 


500 (126) 


550 (126) 


1.3 




556 (125) 


462 (124) 


484 (124) 


547 (126) 


600 (126) 


1.4 




658 (125) 


617 (126) 


525 (124) 


587 (125) 


650 (126) 


1.5 






767 (125) 


635 (125) 


628 (125) 


699 (126) 



Table 2: For different values of the singlet VEV (s — 5 — 10 TeV) corresponding to M z > ^2 — 3.8 TeV, 
the effect of rising the lower limit on the gluino mass between nig = 1 — 1.5 TeV on fine tuning is shown. 
Next to every fine tuning value, the corresponding Higgs mass ( in Ge V) is shown between parentheses. 
The dash means there's no m h <~ 124 — 127 GeV found in the scanned parameter space. 
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